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Abstract: We compare the matrix model and integrable system approaches to calculating the 
exact vacuum structure of general M = 1 deformations of either the basic M = 2 theory or its 
generalization with a massive adjoint hypermultiplet, the M = 2* theory. We show that there 
is a one-to-one correspondence between arbitrary critical points of the Dijkgraaf-Vafa glueball 
superpotential and equilibrium configurations of the associated integrable system. The latter 
being either the periodic Toda chain, for M = 2, or the elliptic Calogero-Moser system, for 
M = 2*. We show in both cases that the glueball superpotential at the crtical point equals 
the associated Hamiltonian. Our discussion includes an analysis of the vacuum structure of the 
M = 1* theory with an arbitrary tree-level superpotential for one of the adjoint chiral fields. 



1. Introduction 



There are a number of ways of investigating the vacuum structure of an Af = 1 supersymmetric 
gauge theory. In this paper we have in mind two such theories, the basic Af = 2 theory and 
also its generalization with a massive adjoint hypermultiplet (known as the Af = 2* theory), 
both with gauge group U(iV), deformed by an arbitrary superpotential for the adjoint chiral 
multiplet. We will focus on two techniques, based on integrable systems [1-4] and on matrix 
models [5-8]. The other closely related approach involves formulating the problem in terms of 
Seiberg-Witten theory [9,10]. 

In the integrable system approach the vacua are determined by the extrema of the conserved 
quantity associated to the Af = 1 deformation. In other words, the vacua correspond to the 
equilibria of the associated flow. On the matrix model side the vacua are determined by 
the Dijkgraaf-Vafa glueball superpotential. In this note we show that there is a one-to-one 
correspondence between these approaches and at an equilibrium point the spectral curve of the 
integrable system is equal to the auxiliary Riemann surface of the matrix model, extending 
the results of [21] for the massive vacua of the Af = 1* theory. For these special vacua, the 
equilibria are stationary with respect to all of the flows of the integrable system signalled by 
the fact that the spectral curve degenerates all the way to genus one. For a general vacuum this 
will not be the case and the spectral curve will only partially degenerate. In order to complete 
the proof we show that at a critical point the value of the glueball superpotential equals the 
Hamiltonian in the integrable system. 

The relation between integrable systems and matrix models has also been investigated 
in [11] for the basic Af = 2 case. As we shall see, our approach is rather different, although 
ultimately must be related. Our approach has some overlap with the recent paper [12] which 
considers, in the context of the Af = 2* theory, the relation between the Riemann surfaces of 
the matrix model and the Seiberg-Witten curve for the Af = 2* theory constructed by Donagi 
and Witten [13]. 



2. The Basic Af = 2 Case 



First we describe the integrable system approach (for a general review of integrable systems 
see [14]). The Coulomb branch of the Af = 2 theory with gauge group U(iV) is identified with 
the moduli space M- m t of the spectral curve E int of the periodic Toda chain, or afline 
Toda, integral system. The spectral curve is the Riemann surface defined by the equation 

F(x,z) = det(xl + L(z)) = , (2.1) 
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where L(z) is the N x N Lax matrix which depends on the canonical variables {pi, qi} via 
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(2.2) 



In the present context, the canonical variables are considered to be complex. The N independent 
flows of the system can be written in Lax form 



d tt L{z) = [M i ,L{z)] , 
where Mj are N x N matrices. It is then clear that quantities of the form 

H = TrW(L(z)) , 



(2.3) 



(2.4) 



for an arbitrary polynomial function W(x), are conserved. In particular, there are N indepen- 
dent conserved quantities that can be taken to be 



Hi = TrL\z) 



i = l,...,N . 



(2.5) 



Note that the {Hi} are independent of z up to a additive constant for i = N. Any quantity 
of the form (|2.4j ) is some function of the {Hi] and will generate a flow d t that is some linear 
combination of the basic flows dt v Note that the basic conserved quantities {Hi} play the role 
of coordinates on .Mint, the moduli space of Ei nt . 



The curve (2.1) can be written more explicitly as 

N 

F(x, z) = Y[(x -di)+z + k 2N z~ l = 

i=l 



(2.6) 



where the 5j = cii(Hj) are an alternative set of coordinates on .Mint- The curve can be written 
in hyperelliptic form by defining 



N 



2z + Y[( 



(2.7) 



i=i 



in terms of which it takes the form 
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i=l 



This has the form of a double-sheeted cover of the complex x-plane with iV cuts joining the 
sheets. Consequently S int is a genus — 1 Riemann surface. 
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The centre-of-mass motion associated to YliPi an d YliQi ls trivial and we can ignore it. 
The remaining N — 1 conjugate angle variables of the integrable system are naturally identified 
with a point ipj, j = 1 . . . , N — 1, in the Jacobian torus of the Riemann surface S; nt . The 
Jacobian torus is defined as follows (for reference on Riemann surfaces see [15]). First we 
choose a canonical set of 1-cycles on £ int , (Aj,Bj), with intersections Aj ■ A k = Bj ■ B k = 0, 
Aj ■ B k = 5jk, j,k — 1, . . . , N — 1. Let uij be the associated set of N — 1 holomorphic 1-forms 
(abelian differentials of the 1st kind) normalized so that § A .oJ k — 5j k . The period matrix of 
Si nt is the N — 1 x N — 1 matrix with elements 

T jk = (h LU k . (2.9) 

Jb 3 

The Jacobian torus consists of points ipj £ C w_1 with the identifications 

4>j ~ ipj + rij + T jk m k , nj,m k £ Z . (2.10) 



We are interested in the flow generated by the conserved quantity of the form H in ( |2.4|) 
where W(x) is some polynomial of degree n. Notice that W(x) itself will depend on z in such 
a way that H is z- independent. For any choice of Hamiltonian H, the associated dynamics is 
linear in the Jacobian. In other words, for each Hamiltonian H, there is a linear flow 

il> j (t)=w j t + ij) j (0) , (2.11) 

where zuj are the angular velocities which just depends on E int . It is determined by the unique 
meromorphic 1-form Q on S int via 

Wj = & Q (2.12) 
Jb 3 

normalized by 

f Vl = (2.13) 

and which is holomorphic on Si nt — P±, where P± are the two points at x = oo on the upper 
and lower sheet. The singularities at P± are specified uniquely by the conserved quantity H in 
the following way. For each polynomial W(x) in ( |2.4j ) there is a unique polynomial W(x) of 
the same order, for which 

^lim Q(P) = ±d(W'(x) + 0(1/ x)) . (2.14) 

Since W(x) is some fixed polynomial and doesn't depend on z by choice, it must be the case that 
W(x) depends on z in such a way that the associated Hamiltonian H in (|2.4j) is ^-independent. 
We will make the relation between W(x) and W(x) more explicit later. 
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We now consider how to use the integrable system to find the vacuum structure of the 
deformed M = 2 theory. In the field theory, the deformation is achieved by adding a tree-level 
superpotential 

-^—TrW(&) (2.15) 

9ym 

to the bare Lagrangian of the M = 2 theory. Here, <3? is the M = 1 chiral superfield of the 
M = 2 theory. In the low-energy effective theory the deformation turns on a potential on 
the Coulomb branch. One can approach the vacuum problem directly in the four dimensional 
theory by using Seiberg-Witten theory. The Seiberg-Witten curve is the spectral curve of the 
integrable system and the Coulomb branch is parameterized by the conserved quantities {Hi}. 
Vacua which survive breaking to M = 1, correspond to special points on the Coulomb branch 
where a number dyons become massless and condense after breaking to M = 1 [9, 10]. However, 
in order to relate the problem directly to the integrable system it turns out to be more useful 
to compactify the theory on a circle to three dimensions [1-3]. This is because the dimension 
of the Coulomb branch is then doubled by the addition of the Wilson lines and dual photons 
of the unbroken U(l) gauge group. The resulting Coulomb branch of the three-dimensional 
theory is identified with the (complexified) total phase space of the integrable system where we 
have not only the Hamiltonians but also the conjugate angle variables. 

It has been shown [1-3] in the context of the M = 2* theory, and its quiver generaliza- 
tions, that the effect of the M = 1 deformation can be captured — including all the quantum 
corrections — in the three-dimensional theory by taking the superpotential on the Coulomb 
branch of the three-dimensional theory to be the conserved quantity H in fl2.4j ). We expect 
these facts to be true in the basic Af = 2 theory as well, since this theory can be obtained by 
taking the large mass limit of the M = 2* theory. This philosophy has also been advocated 
in [11]. 

Putting aside the usual caveats, supersymmetric vacua are critical points of the superpo- 
tential and are therefore points in the complexified phase space which are stationary under the 
flow d t generated by the Hamiltonian H. The value of the superpotential, i.e. H, at the critical 
point is then valid in the four-dimensional limit. So from the integrable system point-of-view, 
as first pointed out in [2], we have to identify equilibrium points of the complexified integrable 
system under the flow dt generated by H. This means that the vacua correspond to points 
in the moduli space where the vector of angular velocities Wj vanishes; in other words, at the 
equilibrium point, due to fl2.12| ) and ( [2.13j ), it follows that 



Q = <h n = . (2.16) 

This implies that there exists a meromorphic function G on E int such that 

Vl = dG (2.17) 
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with singularities at P± = P of the form 



lim G(P) = ±W'{x) + 0(l/x) . (2.18) 
P— *p± 

So to summarize: the supersymmetric vacua correspond to points in the moduli space M.- Ul t 
where there exists a meromorphic function G on Ej nt with a particular pole structure ( |2.18j ) at 
P±. Notice from G we can form the meromorphic function G + W'(x) whose only singularity 
is a pole at P + . Now we can bring the Riemann-Roch Theorem to bear on the question of the 
existence of G. First of all, P + is not a Weierstrass point 1 of S int , since the latter is hyperelliptic 
and the Weierstrass points are located at the 2N branch points. Therefore, the existence of a 
meromorphic function with a pole of order n at P + requires that S int degenerates to a surface 
of genus < n. In other words, in the hyperelliptic form (|2.8|) , it must be that 

N-n 

y 2 = y 2 l[(x-r t ) 2 , (2.19) 

i=l 

where 

2n 

y 2 = X\{x- Sl ) (2.20) 

1=1 

describes a surface of genus < n. In particular, the meromorphic function G is identified with 



y. The pole structure of G (|2.18|) then fixes y uniquely to be 

f = W\xf + f(x) , (2.21) 

where f(x) is some fixed polynomial of degree n — 1. Note that at the equilibrium point where 
S mm degenerates to a surface of lower genus n — 1, the Jacobian must also degenerate. This 
signals the fact that the equilibrium point of dt is also stationary with respect to other flows. 
In general if the reduced curve has genus g then there will be iV — g — 1 stationary flows. 

These facts dovetail completely with the matrix model approach [5] and, for that matter, 
the Seiberg-Witten theory approach [9,10]. Without reviewing the matrix model approach, we 
simply note that for an M = 1 deformation described by W(x) in (|2.15 ), the solution of the 
matrix model involves a hyperelliptic Riemann surface of the form 

f = W '{xf + f(x) , (2.22) 

where f(x) is a polynomial of degree n — 1. This curve is manifestly identical with the spectral 
curve of the integrable system at the equilibrium point under the flow associated to H. 



1 The Weierstrass points {Pi} are the finite set of points on a Riemann surface for which for each Pi there 
exists a non-trivial meromorphic function with a singularity only at Pi with an order less than or equal to the 
genus. 
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To complete the picture, we now prove that the value of the Hamiltonian H that generates 
the stationary flow is stationary is equal to the critical value of the Dijkgraaf-Vafa superpo- 
tential. The fact that we can find the Hamiltonian that generates the flow described by the 
abelian differential Q rests on our knowledge of the conjugate action variables the the angle 
variables [25-27] (the elliptic Calogero-Moser case was first considered in [28]. The conjugate 
angle variables a,j, j = 1 . . . , N — 1, are given by integrals of a certain abelian differential of the 
3rd kind A: 

dj = <p A . (2.23) 



A is precisely the "Seiberg-Witten" 1-form for the U(iV) theory, as deduced in [29,30]. The 
defining property of A is that 

t^-A = 8 jk . (2.24) 
We will take 

A = - \og(y + P)dx , (2.25) 



which is the appropriate form when the derivative in ( |2.24| ) is understood to be taken at fixed 
x. This is necessary in order that the a^-derivative can be pulled out of the a;-integral below. 
Defining 

we see that ttj is an abelian differential of the 3rd kind normalized by 

7r fc = 5jk (2.27) 

and with simple poles at P±. We now apply Riemann's bilinear relation to the abelian differ- 
entials TCj and Q: 



N—l 

-V i SI <b 7Tj — (p SI <b TTj = Hesp + (W'(x)iTj) — Hesp_(W (x)iTj) , (2.28) 

2m k=1 JA k JB k JB k JA k 



where we have used the asymptotic forms for Q at P± in ([2.14 ). Now we use the normalization 



conditions (p. 13 ) and ( [2. 27 ), along with the definition of the angular velocities ( 2.12|) and the 



fact that the contribution from P_ is minus that at P + , to arrive at 

Wj = — Resp, (W'(x)itn) 

47U 



d 

— Res P+ {W\x) \og{y{x) + P{x))dx) ( 2 .29) 

-^-Res P+ ( W(x)^ldx' 
Odj V y{x) J 
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where in the last line we integrated by parts. Now since aj is canonically conjugate to the angle 
ipi, this means that the Hamiltonian which generates the flow associated to Q is 

H = -Res P+ (W(x)^-rdx) . (2.30) 
+ V y(x) J 

It is straightforward to show that if W(x) is a polynomial of order less than N then 

H = TrW(L(z)), (2.31) 

where L(z) is the Lax matrix. So in this case we can identify W(x) and W(x). If W(x) has an 
order > iV then W(x) will be z-dependent. For instance, for W(x) = x N 

W(x) = x N + (-i f (z + A 2N /z) . (2.32) 

What is particularly nice about the result ( |2.3U| ) is that at a critical point it agrees precisely 
with the critical value of the Dijkgraaf-Vafa superpotential. In order to see that, one simply 
deforms the contour around infinity to a sum of contours around the cuts on the degenerated 
surface. We then identify the differential P' dx/y with the resolvent of the field $ defined 
in [31,32] 

T(x) = Tr^-. (2.33) 
x — <P 

The critical value of the superpotential is then 



W(x)T(x) (2.34) 

' cuts 

which is equal to ( |2.30| ). Hence we find perfect argument between the two distinct approaches 
for calculating the vacuum structure. 



3. The M = 2* Case 



We now apply the same philosophy established in the basic M = 2 case to M = 1 deformations 
of the Af = 2* theory. 

The Coulomb branch of the U(iV) Af = 2* theory is identified with the moduli space of the 
spectral curve S int of the iV-body elliptic Calogero-Moser integrable system [13,16]: 

F(v,z) = det(vl + L(z)) = , (3.1) 

where the N x N Lax matrix L(z) has components 

L t3 (z) = Pl 5 t3 + m(l - 5 l3 ) a \\-? l + q >\ e^-^ . (3.2) 

a(z)a(q i -q j ) 
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We denote the iV-dimensional moduli space of the curve by .Mint- The quantity v is a mero- 
morphic function on S int with iV simple poles at the pre-images of z = with residues 

m(JV-l,-l,-l,...,-l) • (3.3) 

Notice that one of the points, which we denote Po, is distinguished by the fact that the residue 
is m(N — 1). In the Type IIA brane construction of Witten [17] Po is the position of the 
NS5-brane. From a field theory perspective, m is the mass of the adjoint hypermultiplet and r 
is the bare complexified gauge coupling. 



F(v, z) is a polynomial of degree N in v whose coefficients are elliptic functions on the 
torus E T with complex structure r: 

N 

F(v,z)=Y,M*)v i > (3-4) 

where 

fi(z + 2m) = fi{z + 2mr) = f t (z) . (3.5) 

The spectral curve S int describes an A^-sheeted cover of the torus E T joined by branch cuts to 
make a higher genus surface. The number of branch cuts corresponds to the number of zeros of 
d v F(v, z). Since the latter is a meromorphic function on Ej nt the degree of its zeros is equal to 
the degree of its poles. It follows from (|3.3j ) that d v F(v, z) has A^ — 1 simple poles and a pole 
of order AT — 1 at the A" pre-images of z = 0. Hence, there are 2(A" — 1) branch cuts and the 
Riemann-Hurwitz Theorem gives the genus of S int as N. We can view the surface as A" copies 
of the torus E T glued together by A^ — 1 tubes to make a genus A^ surface. 

Any quantity of the form 

H = TtW(L(z)), (3.6) 

for a polynomial function W(x), will be conserved. A basis for the Hamiltonians, and so a set 
of coordinates for A^ int , is provided by 

Hi = Trl?(z) , i = l,...,N, (3.7) 

where z takes a fixed value ^ 0. The A^ conjugate angle variables are, as before, associated 
with a point in Jacobian of S int . The dynamics is linear in the Jacobian [14, 18, 19] and for 
the flow corresponding to an arbitrary Hamiltonian H in (|3.6|) , we have 

^■00=^ + ^(0), (3.8) 

where Wj is the vector of angular velocities associated to H. This quantity is determined by 
the unique meromorphic 1-form Q on Ei nt with 

n = o , Wj= I n (3.9) 
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such that it is holomorphic on T,- mi — P$ with a given singularity at Pq determined by H. Let x 
be a coordinate in the neighbourhood of Pq with x(Pq) = oo, then the singularity of Q has the 
form 

\im Q(P) = d(U(x) + 0(1/2)) , (3.10) 

where U(x) is a polynomial in x of the same degree as W(x). In our philosophy, U(x) will be 
fixed uniquely by the M = 1 deformation, as we shall see later. This then fixes the Hamiltonian 
H and hence the function W(x). In particular, since U(x) is by choice z- independent, W(x) 
must be ^-dependent in such a way that H in (|3~6| ) is independent of z. 



Following the logic of the last section, we now use the structure of the integrable system 
to find the vacuum structure of the M = 2* theory deformed to M = 1*. The deformation is 
achieved by adding a tree-level superpotential 

-^—TrW(&) (3.11) 

9ym 

to the bare Lagrangian of the M = 2* theory. Here $ is massless adjoint chiral superfield of the 
M = 2* theory. As before we compactify to three dimensions and identify the Coulomb branch 
of the theory with the (complexified) total phase space of the integrable system. 

The deformation (|3.11| ) is captured exactly by the superpotential on the Coulomb branch of 
three-dimensional theory which is equal to the Hamiltonian H whose flow was described above. 
Note that in general the two polynomial functions W(x) and W(x) are not equal, however, 
they are of the same order and are uniquely related to one another — at least implicitly — as we 
shall see later. Supersymmetric vacua are critical points of the superpotential and are therefore 
points in the complexified phase space which are stationary under the flow d t generated by the 
Hamiltonian H. Consequently, at these point the associated angular velocities must vanish; 
hence, 



n = f n = . (3.12) 

Aj Jbj 

This implies that there exists a meromorphic function G on S int such that 

Vl = dG, (3.13) 
with, from (|3.10 ), a singularity at Pq, of the form 



lim G(P) = U(x) + 0(l/x) . (3.14) 

P— >Po 



If W(x) is a polynomial of degree n then so are both W(x) and U(x). Generically, according 
to the Riemann-Roch Theorem, for G to exist, S int must degenerate to a surface of genus n — 1. 
However, if P happens to be a Weierstrass point of S int then there may exist additional vacua 
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where Si nt has genus > n — 1. In particular, n would have to be in the "non-gap" sequence at 
Pq. Unlike in the hyperelliptic situation described in the last section these exceptional cases do 
indeed occur. An example occurs in the U(3) theory with the simplest quadratic deformation 
W(x) = x 2 . In this case, there is a vacuum described in [13] where the surface degenerates from 
genus 3 to 2 but not all the way to genus 1 as would be required if P was a generic point on 
the surface. In particular, for the quadratic deformation, as was first pointed out in [2] all the 
vacua must be described by a degeneration to a hyperelliptic surface since it is only on these 
surfaces that there exists a meromorphic function with a double pole. In particular, P , being 
a Weierstrass point, must lie at one of the branch points in the two-sheeted representation. We 
will leave a more in-depth discussion of these exceptional cases to future work. 

Before we move on to describe the matrix model approach to the same problem, we note 
that we can map the curve S int into the complex plane parameterized by x via 

x = iv(z) + im(c(z) - • (3.15) 

Here, ((z) is the Weierstrass ^-function which is a quasi-periodic function on E T : 

C(z + 2m) = C(z) + 2((m) , ((z + 2mr) = ((z) + 2((mr) . (3.16) 

In the rr-plane, P is mapped to x = oo. Note that x is multi-valued on E int because, although it 
is periodic around the pre-images Oj of the a-cycle of the base torus E T , it picks up an additive 
constant im around the pre-images hi of the 6-cycle of the base torus. So restricting x to a single 
sheet there are N pairs of cuts C~ = [ki, k] and = [fcj + im, ^ + im] which are identified. In 
this picture the genus iV surface S int is realized as the complex x-plane with N handles formed 
by identifying the top/bottom Cf with the bottom/top of C~ . In general a Riemann surface of 
this form has 2N complex moduli which we can take to be the positions of the ends of the lower 
cuts {ki,li}. We denote the moduli space of these surfaces as M.. Clearly the moduli space 
of TV-fold covers of the base torus E T , Aimt, is only a complex iV-dimensional subspace of this 
larger moduli space. Note that when the surface S int degenerates, cuts in the x-plane merge. 
Note that two pairs of cuts can annihilate in two distinct ways. Either the bottom cut of one 
pair merges with the bottom cut of another pair to leave a single pair of the same kind, or the 
top cut of one pair merges with the bottom cut of another pair resulting in another pair of cuts 
now separated by 2im rather than im. By iterating this procedure we see that degenerations 
of S int are manifested in the complex x-plane by pairs of cuts which join to form handles which 
can be separated by any integer multiple of im. 

Before we discuss the matrix model approach let us consider the vacuum structure from 
the point-of-view of the tree-level superpotential: 

W = ^Tr(i$[$+, $-] + m$ + $- + W($)) . (3.17) 
9ym 
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The case with a quadratic superpotential W($) was considered originally in [13], here, we 
present the generalization for arbitrary polynomial functions W(§). At tree-level, the super- 
symmetric vacua can be found by solving the F-flatness conditions 

[$, $+] = i m $+ , [$, $-] = -im<$>- , [$+, $-] = zW'($) (3.18) 

modulo complex gauge transformations. Using the symmetries we can diagonalize $ and the 
build up solutions from a series of irreducible blocks. In such a block of size p {$ ± , $} have 
the same non-zero elements as { J^, J3} of the irreducible representation of SU(2) of size p. In 
particular 

$ = zmJ 3 + Al. (3.19) 
The parameter A is then determined by demanding 

p 

Tr block [$ + ,$1 = = ^IU / (A + ^(p+l-2j)) • (3.20) 

If W(x) is a polynomial of degree n, there are n — 1 possibilities for A. Hence, a general vacuum 
corresponds to the data 

9 

{rij, pj, Xj, j = 1, . . . , g} , N = ^2 n jPj (Pi ^ Pj when x i = x i) > ( 3 - 21 ) 

3=1 



where rij is the number of blocks of size Pj associated to one of the n — 1 roots Xj of ( |3.20j ). 
The unbroken gauge group in this vacuum is 

9 

n u («i) ( 3 - 22 ) 

and in particular the number of abelian factors is g. Quantum mechanically we expect in the 
infra-red that the non-abelian parts of the gauge group confine to leave an abelian theory with 
gauge group U(l) 9 . Note that the maximal value of g is N, obtained when the potential W(x) 
is a polynomial of degree > iV and Pj — rij — 1, j — 1, . . . , N. In this case the unbroken gauge 
group is U(l) JV . This is the M = 2* analogue of the vacuum considered in [10] which can track 
the Coulomb branch of the M = 2* theory and be used to extract the Donagi-Witten curve 
from the matrix model (as recently considered in [12]). The minimal value of g = 1 is obtained 
when pifii = N. In other words when p± is an integer divisor of N. These are the "massive" 
vacua considered in [1,20,21] (note that in the U(iV), as opposed to the SU(iV) theory, there 
is an unbroken U(l) factor and so strictly speaking the vacua are not massive). 

Now we turn to the matrix model approach to calculating the exact superpotential. Ac- 
cording to Dijkgraaf and Vafa we consider a matrix model whose matrix variables are associated 
to the chiral superfields of the theory and whose action is the F-term. In the M = 2* — > 1* 
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theory this yields a matrix model with a partition function (see [7, 20-24] for previous work on 
matrix models and the Af = 1* theory and its generalizations): 

Z = y[d$ + ][d$-][d$]exp-^J lr Ir («$[$+,$-] +m$+$- + W($)) . (3.23) 

Since $ ± appear Gaussian we may integrate them out: 

r e -gs 1 TrV(<s>) 
Z -./^ de t (Ad Jt + lm) - < 3 ' 24 > 

In order to avoid confusion, we will suppose that the matrices have size N. In order to implement 
the Dijkgraaf-Vafa procedure to the vacua described above, we need to solve the matrix model 



Eq. ( 3.24 ) in the large iV-limit around the saddle-point corresponding to the classical solution 
for where <3> takes its tree-level form with n, replaced by arbitrary variables rij which individually 
tend to infinity. As usual in the large-n^ limit the eigenvalues of $ spread out from their classical 
values along cuts on the complex eigenvalue x-plane. In other words, for each j = 1, . . . , g there 
will a set of pj cuts which form a group, each element of which being the image of the lower 
one under translations ikm, k — 1, . . . ,pj — 1. Each group is located in the vicinity of Xj. So 
for each j = 1, . . .,g there are 2 complex parameters which one think of as the two ends of 
the lower cut. The density of eigenvalues p(0) only has non-zero support along the cuts in the 
x-plane. Moreover, for a each j = 1, . . . , g the density of eigenvalues along each of the pj cuts in 
the group is equal. The saddle-point equation in the large- iV limit is most conveniently written 
in terms of the resolvent function 

w(x) = f ^rd(f) , ! p{<f))d<j) = 1 . (3.25) 

J cuts V J cuts 

The resolvent u(x) is an analytic function on the complex x-plane whose only singularities are 
branch cuts along the cuts where the eigenvalues are located for which the discontinuity across 
the cut gives the eigenvalue density 

u{4> + * e ) ~ ^(0 — ^ e ) — — 27rip(0) , G cuts . (3.26) 

In this, and following equations, e is a suitable infinitesimal regulator. The saddle-point equa- 
tion expresses the condition of zero force on a test eigenvalue in the presence of the large-TV 
distribution of eigenvalues along cuts: 

W'(d>) 

— - — = u)((f) + ie) + u)((f) — ie) — u{<p + im) — cu((fi — im) , G cuts , (3.27) 

where S = g s N is the 't Hooft coupling. This equation can be re-written in terms of the useful 
function 

G(x) = U{x) + iS{u{x + if) - u{x - , (3.28) 
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where U(x) is a polynomial in x of the same degree as W(x) such that 

W'{x) = -iU(x + if) + iU(x - if) . (3.29) 

It turns out that G(x) has a somewhat simpler analytic structure than u(x). For a given group 
of pj cuts the situation is described in [21]. In taking the difference of the resolvents in ( p.28|) 
most of the cuts cancel to leave only a pair of cuts associated to each group [kj,lj] and [kj + 
impj, lj + impj]. Moreover the saddle-point equations simply imply a gluing condition between 
each pair of cuts where the bottom/top of the lower cut is identified with the top/bottom of 
the upper cut. This naturally defines a Riemann surface of genus g since each pair of cuts 
once glued corresponds to one handle. In fact it should not escaped the reader's notice that 
this Riemann surface is precisely an example of the family of surfaces M. defined earlier in the 
context of the integrable system. 

However there is more structure since the Riemann surface has to admit a meromorphic 
function G with a singularity at the point P Q , x(P ) = oo, with the following pole structure 

lim G(P) = U(x) + 0(l/x) . (3.30) 

P— >Pq 

So the saddle-point equation of the matrix model boils down to the existence of a genus g 
Riemann surface S mm in the space M. which admits a meromorphic function G whose only 
singularity is at the point Pq with the specific pole structure (|3.30|) . 

The question is how many moduli does the surface S mm have? Potentially there are 2g 
complex moduli as pointed out above. However, there are additional constraints arising from 
the requirement that the function G exists on S mm . Generically, by the Riemann-Roch theorem, 
such a meromorphic function only exist on a surface of genus g < n, in which case there would be 
n — g non-trivial meromorphic functions with a pole of this order or less. Therefore prescribing 
the asymptotic form (|3.30|) amounts to g non-trivial conditions of the surface S mm , leaving a g 
complex dimensional subspace of moduli. The g moduli of these surfaces can be described by 
the g quantities 

s i = -^-i G ( x ) dx > ( 3 - 31 ) 

where Aj is a contour which encircles the lower cut [hi, lj] of each pair. Notice that the number 
of matrix model eigenvalues associated to each j = 1, . . . , g is equal to pjSj/g s and furthermore 

9 

Sj = g s hj , S = g s N = ^PjSj . (3.32) 

Each of the moduli Sj will become a field of the Dijkgraaf-Vafa superpotential identified with 
the glueball superfield of the unbroken U(n 3 -) factor of the gauge group. 
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The second ingredient required to construct the Dijkgraaf-Vafa superpotential is the vari- 
ation of the genus zero free energy JF with respect to Sj. Following [21], this is equal to 

< ^' < —i (f G(x)dx , (3.33) 



9S j jBj 

where Bj is the conjugate cycle to A,-, in other words, goes from the lower cut to the upper cut 
of a pair. The quantity ( ^.33[ ) can be interpreted physically as the variation of the genus zero 
free energy of the matrix model in transporting pj eigenvalues in from infinity and placing one 
on each of the pj cuts in a group (so as to maintain the same density along each of the cuts in 
the group). We conjecture that the generalization of the Dijkgraaf-Vafa glueball superpotential 
for the vacuum in question is then 

dF 



W eS (Sj) = £ fadi ~ ZniPtrS,) , (3.34) 
j=i J 

where now the rij are the physical rather than matrix model quantities. 



A critical point of W c s corresponds to 



£%J^ = 2 ^*- (3 - 35) 
j=i 

This equation can be written in a more suggestive way by noticing that 

1 d 

ujj = - — ——G(x)dx , (3.36) 

are a basis for the g holomorphic 1-forms on E mm since the singular part of G(x)dx at x = oo 
is manifestly independent of the moduli {Sj}. Furthermore, the iOi are normalized so that 



Hence 



u k = 5 jk ■ (3.37) 



2ixi (p Uk = 2'niTjk , (3.38) 



dSkdSj J B . 

where Tjk is the period matrix of S mm . Consequently the critical point equations are 

9 

^2rijTj k = rp k . (3.39) 
j=i 

These equations are precisely the conditions that S mm is an iV-sheeted covering of the base 
torus E T . In order to prove this we need to find a map from S mm to E T which covers the latter 
N times. For P 6 S mm , the map is simply 

pP 9 

z(P) = 2ni / 2_, n j UJ j m od 27ri,27rrr , (3.40) 

' 3=1 
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where P' is a fixed but arbitrarily chosen point of S mm . Since A" = J2j n jPj this map covers 
E T N times. In particular, z is identified with the same quantity on the integrable system side. 
Similar ideas have been expressed in [12] for the case when g = N, rij — p 3 — 1 corresponding 
to the vacuum with unbroken U(l) symmetry. 

Now we compare the integrable system with the matrix model. The first point is that 
the moduli space .Mmt and the union of the moduli spaces of the matrix model for different 
vacua -M mm are both subspaces of the moduli space M. defined earlier. The subspace .Mint 
contains the surfaces which are iV-fold branched coverings of the torus E T and the subspace 
Ai mm contains surfaces which admit a meromorphic function G with a fixed pole structure at 
Pq determined by the Af = 1 deformation. The critical-point condition in the integrable system 
approach boils down to the existence of the meromorphic function G while on the matrix model 
side it boils down to the constraint of being an iV-fold cover of the torus E T . In other words, 
the vacua correspond to points of intersection between M.^ and the different components in 
-M mm . In this picture we must identify the polynomial U(x) on both sides of the story which 
then implicitly determines the relation between W(x) and W(x). Notice for the case of the 
vacuum with g = N and unbroken U(l) N we can, generalizing the situation in the basic M = 2 
theory [10], extract the Donagi-Witten curve from the matrix model model as was pointed out 
in [12]. 

Now we turn to the question of the value of the superpotential at the critical points cal- 
culated using the two methods. We start by finding a more explicit expression for the critical 
value of the superpotential on the matrix model side by the following manipulations: 

9 



W e fi = —i (rij ® G{x) dx — rpj (b G{x) dxj 

j=l J B i J a j 

rij ( (b G(x)dx — Tjk (t G{x) dx 

7 = 1 fe=l 



3=1 
9 

-l 

3~- 

-2tt Resp (U(x)z dx) . 



(3.4i; 



In the above, to reach the second line we used the critical-point equations ( |3.39| ) and to reach 
the final line we applied a Riemann bilinear relation, used the fact that dz = Ylj n j UJ j an d that 
G(x) can be replaced with U(x) in the vicinity of Pq. A similar expression was derived in [12]. 

Now from the integrable system side. Once again we use the fact that action variables 
conjugate to the angle are integrals of the Seiberg-Witten differential: 

aj = i X , X = vdz . (3.42) 

Since the residue of A at P is independent of the moduli a i; the derivative of A with respect 
to di is the holomorphic 1-form u>i. However, with the Seiberg-Witten differential in the form 
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v dz we must be careful to differentiate at fixed z otherwise singularities arise: 

It turns out that this is not a conventient form. One can eaily verify from the fact that 
F(v, z; di) = and from the relation between x and v and z, that 

„,= (*) dx. (3.44) 

V Odj J z \ Odj J v V OCL; / i 



x 



The final expression here is the one which will be most convenient. Applying a Riemann bilinear 
relation to Q and u>i one arrives at the following expression for the angular velocities 

d 

Wj = — Resp (U(x)z dx) . (3.45) 

Hence the Hamiltonian which degenerates the flow described by the abelian differential Q is 

H = -2vr Res Po (U{x)z dx) (3.46) 

which is in perfect agreement with the matrix model result fl3.41| ). In a sequel we shall show 
how to write ( |3.46| ) in terms of the Lax matrix of the integrable system and hence has a function 
of the positions and momenta. 

Obviously important questions remain, the most interesting concerning the exceptional 
vacua for which the surface has a genus greater than or equal to the order of W(x) and for 
which Pq must be a Weierstrass point. 



I would like to thank Harry Braden, Jan de Boer, Nick Dorey, Prem Kumar and Kazutoshi 
Ohta for discussions. 
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